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They say a picture is worth a thousand words. It is often the case in mathe¬ 
matics that a visualization of the objects we work with is worth a thousand equations; 
visualizations allow us to observe the connections hidden behind layers of rigor and 
equations. In the case of Figure 1, it is worth 67,108,856 equations, all of which are of 
the form p{z) = 0. Here, all p{z) are Littlewood polynomials; Figure 1 is simply the 
plot of polynomial roots. Specifically, Littlewood polynomials are polynomials with 
coefficients of either 1 or -1, and were first studied in depth by J. E. Littlewood during 
the 1950’s. Thanks to recent advances in computing, it is now possible to compute the 
zeros of this class of polynomials, under some restrictions. Surprisingly, when the ze¬ 
ros for all Littlewood polynomials, up to a certain degree, are plotted, (see Figure 1 and 
Figure 2), the zero sets are locally fractal-like and exhibit high degrees of symmetry. 

Motivated by the complex symmetry of the zero set of Littlewood polynomials, we 
examine the symmetry of the set of Littlewood polynomials themselves. In Figure 2 we 
see that the distribution of the zero sets are locally fractal-like. Thus, the complexity 
of the set grows exponentially on a local scale. The computational complexity of such 
sets is to be expected, since computing this set means solving the root finding problem 
for potentially millions of polynomials. We use the relationship between the symmetry 
in the zero set of Littlewood polynomials and in the symmetry in their coefficients to 
reduce the computational complexity of the zero set and to give an alternate character¬ 
ization of the set. Namely, we will study connections between the group structure of 
the coefficients of Littlewood polynomials, the distributions of the zeros and a fractal 
set called a generalized dragon set which will be defined along the way. 

Before we jump into exploring the symmetries of fractals and fractal-like objects, 
we should be clear on what a fractal is. While a universally accepted definition of a 
fractal set has yet to be established in the literature, in the famous words of US Supreme 
Court Justice, Potter Stewart, ”1 know it when I see it”. This very unmathematical and 
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Figure 1: Roots of Littlewood polynomials up to degree 24 in [—2, 2] x [—1.5,1.5] C C. 
Coloration denotes density. 

Image credit: S. Derbyshire 



Figure 2: Close up of Figure 1, first quadrant of C. 

Image credit: S. Derbyshire 


fuzzy manner of identifying fractals and fractal-like objects is perhaps due to fractal 
geometry being a very new branch of mathematics. Fortunately, if one is a bit more 
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selective with the type of fractals one wishes to study, it is possible to provide a solid 
definition. For our purposes, we shall concern ourselves with iterated function systems, 
usually written IFS. An IFS is a finite set of contraction mappings on a complete metric 
space (for this paper the reader may just assume that by complete metric space we mean 
or C). The way one constructs a fractal from an IFS is by iterated applications of 
contraction mappings to a compact set. If we let {/j \ X X \ j = 0,1, 2,..., d}, be 
a set of contraction mappings on a compact space endowed with the Hausdorff metric, 
X, (here, the reader may assume X is roughly or C along with a notion of distance 
between subsets) then, by the Banach fixed-point theorem [7], there exists an attractor 
set 

S = lim Sn^ 


where 


d 

Sn+l = U fASn), 

i=i 


which will be called a fractal. 

In order to solidify the notion of an IFS, we will show the construction of a well known 
fractal set, the Sierpinski Triangle. We let A = C, fi{x) = f 2 {x) = + ^, 

and f 3 {x) = -h (| + be our contraction mappings, and Sq be an equilateral 
triangle region in C with vertices at 0,1, and ^ Then, as outlined above, the set 

Si = fi{So) U f 2 {So) U f 2 {So). Thus, Si is the union of three maps of Sq. The plots 
below visually show this process for up to n = 3. 





In the first plot we see the region Sq. In the second plot, we see Sq U Si with Si 
being a darker shade. Similarly, the darker regions of the third and fourth plot are S 2 
and S 3 respectively. In short, what this process is doing is taking 5'n-i, scales it down, 
makes three copies, and then applies different translations to each copy to produce Sn- 


Another excellent, and much more relevant example of an IFS is the Littlewood 
polynomials themselves! Recall that the Littlewood polynomials are polynomials with 
coefficients of ±1. With that in mind, let us see how these polynomials may be con¬ 
structed through function iteration. Let /+ (x) = 1 + zx and /_ (x) = 1 — zx be 
contraction mappings from C to C. If we start to iterate these two contractions we get 
something like 

/+(/-(/+(/+(0)))) = l-z-z‘^ +z^. 


Notice that no matter how we compose /_ and /+, we shall always get a polynomial 
with coefficients of ±1. However, the set of Littlewood polynomials of some fixed 
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degree d, which we shall denote as 



e{-i,i} 


includes polynomials with —1 as a constant term. Therefore, by iterating /+ and /_ 
we only get half of the Littlewood polynomials, namely, the ones with a constant term 
of 1. Fortunately, since we are concerned with their roots, and the IFS construction of 
Cd is only off by a factor of —1, this IFS construction will suffice. In our first example 
of an IFS, our final object that we ended up with was some geometric set in M?, but 
in this second example, we ended up with a set of functions. This is perfectly okay! 
It is exactly this fractal of functions that can then give us a ’’geometric” fractal. The 
type of fractal object that one gets from this collection of functions is what is known 
as a generalized dragon set. Therefore, a generalized dragon set will be the set of all 
polynomials in Cd with a constant term of 1, evaluated at some z G C. In Figure 3a 
and Figure 3b below, we illustrate this fact below with a plot of 


-^13(^0) = {^(^o) I p e .C13}, 


for fixed 2^0 G C and of its corresponding Generalized Dragon Set. 
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Figure 3a All Littlewood polynomials of degree 13 evaluated 2^0 = 0.48 + 0.45i 
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Figure 3b 14 iterations of generalized dragon set with parameter zq = 0.48 + 0.45i 

Now that we have seen the connection between Littlewood polynomials and gener¬ 
alized dragon sets, we will continue our study of the roots of Littlewood polynomials 
via iterated function systems. Not only does this connection outline where the locally 
fractal-like structure of the roots comes from but it also allows us to exploit the high 
degree of symmetry in the set of polynomials to approximate the roots, which we shall 
denote as 

Dd = {ze<C\ p{z) = 0, Vpe £d} ■ 

A visualization of such a set can be found in Figure 1, as it is D 24 . Naturally, we 
shall use groups, to capture this symmetry. Thus, we wish to find an appropriate group 
structure on 

Two obvious choices for the binary operation on Cd would be either the ordinary 
polynomial product, or polynomial addition. However, we immediatly see that neither 
of those would preserve the requirement that the coefficients be ±1. Instead, we use a 
matrix product from Linear Algebra, the Hadamard product, to induce a group struc¬ 
ture to the set C. 

The Hadamard product, which we shall denote as o, is defined to be component-wise 
multiplication between two matrices of the same size. 

For example, the Hadamard product of two row vectors is given as (ai, a 2 ) o ( 61 , 62 ) = 
(ai * 61 , a 2 * 62 ), where * denotes regular multiplication of scalars. 

Equipping Cd with o gives an Abelian group (that is, a group where the operation is 
commutative) with 2^+^ elements. The closure, associativity, and commutativity of o 
all come from the closure, associativity, and commutativity of * over the set { — 1 , 1 }. 
Furthermore, the identity element will be the polynomial of degree d with all coef¬ 
ficients being 1, and every Littlewood polynomial is its own inverse in this group. 
The reader may have noticed by now that the group (Cd^ o) is naturally isomorphic to 
Z 2 , where 0 denotes the direct sum of groups, and Z 2 is the group of integers 
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mod 2 under multiplication mod 2. Essentially, we pick up a copy of Z 2 for every 
coefficient in a Littlewood polynomial. 

It is known that ^2 is finitely generated, and so Cd must also be finitely 

generated by some subset Sd C with d+l elements. For our purpose, we shall only 
consider Sd to be the set of Littlewood polynomials with only one negative coefficient. 
That is. 
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Now that a group structure on Cd has been established, the existence of a group 
generating set will allow for a coordinate system to be established on Cd{z) = {p{z) : 
p G Cd}. Ultimately, the generating set for Cd will allow us to factor any Littlewood 
polynomial over the Hadamard product into a product of polynomials in the generating 
set. Note that this will not be factoring in the usual sense, i.e. over the canonical 
polynomial product (that is, multiplication of polynomials in the usual sense using 
the distributive property). To get a sense of polynomial multiplication in the sense of 
Hadamard, let us compare it with the usual product of polynomials. For example, let 
p{x) = —x‘^ + X + 1, and q{x) = —x‘^ — x + 1. Then, the Hadamard product of p{x) 
and q{x) will be: 

p{x) o q[x) = x^ — X ^ 1. 

Compare this with the usual way of multiplying polynomials: 

p{x)q{x) = + 1 . 

It will turn out to be the case that by factoring any polynomial in Cd in the sense of 
Hadamard, we shall be able to uniquely write any point in £^(^ 0 ) as a binary linear 
combination of a much smaller subset of Cdiz^). That is, we shall be able to compute 
the evaluation of a large number of polynomials without actually evaluating too many 
of them. 

The following image highlights the points in £i3(0.48H-0.45i) which are the images of 
a polynomial in 5 'i 3 in red and the point which is the image of 613 in green. In short, we 
will show that every point in blue can be expressed as a binary sum of points in red, plus 
a correcting term which will be a scalar multiple of the point in green. Furthermore, 
we will show how if a scalar multiple of the point in green can be expressed as a 
linear combination over { 0 , 1 } of the points in red, then the parameter z used in the 
construction of this image is a root of some Littlewood polynomial. 
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Figure 4 £ 13 ( 0.48 + 0.45i) with 5'i3(0.48 + 0.45i) in green, and 613 (0.48 + 0.45i) 

in red. 


Factoring Function: 

Since the group Cd is generated by some subset, it is only natural to ask how do Little- 
wood polynomials factor over this generating set, and how does polynomial evaluation 
interact with the process of factoring in the sense of Hadamard. In order to gain some 
insight, we need to define a factoring function. Let p e Cd- Then, define a factor 
function, cr(p), to be a subset of the generating set, Sd, of jCd such that the product 
of its elements is p. In short, since Cd is generated by some subset, the factoring 
function a tells us how to write any polynomial in Cd as the Hadamard product of Lit- 
tlewood polynomials in the generating set. Furthermore, we define a measure function 
a \ Cd ^ ^ on Cd such that ^{p) is the number of negative coefficients of p. That is, 
iy{p) = |cr(p)|, where \cr{p)\ is the number of elements in the set cr(p). 

For an example of a and z/, consider the Littlewood polynomial — x‘^ — x + l G £ 3 . 
Then, a(x^—x‘^—x+l) = {x^ —x‘^-\-x-\-l, x+1}, and u{x^—x‘^—x-\-l) = 2. 

The following lemma and theorem give a way of taking the generating set of Cd 
and using it to write any point in £^^(^ 0 ) as a linear combination of points in S'^( 2 ;o). 
Essentially, they will summarize how polynomial evaluation and the Hadamard product 
of polynomials interact with one another. 

Lemma 1 Let Sd be the generating set of Cd with d> 2 and z e C. Then, 

e{z) = ^ ^ g{z). 
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Proof. We begin by observing that e{z) = • Furthermore, 


Thus, 


d+l d 

g^Sd i=i«=o 


d 


{d-i)Y,z\ 

i=0 


e{z) = 


1 

(i — 1 


I] 9{z). 

geSd 


□ 


As an example of Lemma 1, we let (i = 3 and zq = 0.5 + 0.5i. Then by Lemma 1, 

^(^o) = 2 ((~^ 0 +^ 0+^0 + l) + (^0 ~^ 0+^0 + l) + (^ 0+^0 ~^0 + l) + (^ 0 +^ 0 +^ 0 ~l)) 

e(0.5 + 0.5z) = i(1.75 + 0.75i + 1.25 + 0.25i + 0.25 + 0.25z - 0.75 + 1.25z) 
e(0.5 + 0.5i) = 1.25 + 1.25i. 

In short. Lemma 1 gives a way of evaluating the identity element in Cd in an indirect 
way. 

The next theorem then allows us to do the same, evaluate arbitrary polynomials in the 
group jCd in an indirect manner. 

Theorem 2 Letp G £. If o is the Hadamard product then, for all z G C, 
p{z) = {1 - I^{p})e{z) + y] g{z). 

gecr{p) 


Proof Let p G C such thatp = op 2 o * * * op/c where gi G Sd and k = iy{p). Also, let 
and denote the component of p and g respectively. Then, if p^ = 1 it follows 
that gj = 1 for all 1 < i < /c. Thus, 



{k — 1) = k — {k — 1) = 1 = p^. 


Now if = —1, there exists only one gi such that gj = — 1. Thus, 


I - {k-1) = {k-2) - (k-l) = -1 =pL 


j=i 


Therefore, in both cases. 


p^ 



{k - 1). 
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Thus, we conclude that 




\i=i / 

Thus, since the polynomials p, g and e, are all in C[x], it follows that for the polynomial 
functions, one has 

p{z) = (I - v{p))e{z) + g{z). 

gecT{p) 

□ 

As an example of Theorem 2, we let p e Cs such that p{x) = — x‘^ — x -\-l, and 

zq = 0.5 + 0.5i. Then 


p{zo) = (1 - 2)e(zo) + {zq - zl + zo + l) + ( 2:0 + zl- zq + I) 


p(0.5 + 0.5z) = -1.25 - 1.25i + 1.25 + 0.25i + 0.25 + 0.25z 
p(0.5 + 0.5i) = 0.25 - 0.75i. 

It has now been shown how the group structure of Cd, the group of Littlewood polyno¬ 
mials of degree d under the Hadamard product, can be used to establish an addressing 
system on Generalized Dragon Sets, which we have already shown to be what Dd looks 
like on a local scale. That is, any point in the Generalized Dragon set, Cd{z), can be 
expressed as a linear combination of the elements of the generating set, S{z) and the 
identity, e(z), of Cd- This addressing system, taken together with the fact that Dd lo¬ 
cally looks like Cd{zQ) for some zq G C, gives us a sort of local coordinate system on 
Dd. 

This then allows us to essentially slice up Cd{zo) into smaller, more manageable parti¬ 
tions. In mathematics, mathematical objects with an underlying set, such as our groups 
of Littlewood polynomials, are partitioned by using equivalence classes. Formally, 
we begin by observing that 1 < u{p) < d, with u{p) G N. Let Nd{n) = {p G 
• ^(p) = be an equivalence class given by the partition induced by z/. Then, 

\Nd{n)\ = Cfp)), and ULo = U 

To provide a visualization of such a partition on a dragon set, consider the plots below, 
which are the plots of the evaluations of polynomials in Ni 4 {n) at z = 0.48 -j- 0.45i, 
for 1 < n < d. Evaluations of polynomials in Si 4 {z) are shown in red. 
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(a) Evaluation at 0.48+0.45i of Littlewoodpoly- (b) Evaluation at 0A8+0A5i of Littlewoodpoly¬ 
nomials of degree 14 with iy{p) = 2. nomials of degree 14 with iy{p) = 4. 

Figure 3 



(a) Evaluation at 0A8+0A5i of Littlewood poly- (b) Evaluation at 0A8+0A5i of Littlewood poly¬ 
nomials of degree 14 with o^p) — 6. nomials of degree 14 with v^p) — 7. 

Figure 4 
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(a) Evaluation at 0.48+0.45i of Littlewoodpoly- (b) Evaluation at 0.48+0.45i of Littlewoodpoly¬ 
nomials of degree 14 with iy{p) = 11. nomials of degree 14 with ^{p) = 14. 

Figure 5 


From these plots, we observe that Cd{z) can be approximated by evaluating poly¬ 
nomials in Nd{k) C Cd, where k is some fixed constant in N. When k is chosen to be 
small, the number of polynomial evaluations is large and when k is chosen to be closer 
to d the number of polynomial evaluation is small. We observe in the above plots that 
the evaluations of polynomials in Nd{k) is a. refiection of the evaluation of polynomi¬ 
als with Nd{d — k). Thus, an approximation of this type will require at most 
polynomial evaluations. Indeed, we observe in the above plots, for d = 14, choosing 
k = 7 provides the best approximation. Therefore, we observe that even the “largest” 
approximation of this type will be a set with (pd-|) elements as opposed to Cd{z) which 

will have 2^+^ elements. For a comparison, the following plot shows the difference in 
cardinalities between Cd{zo) and the best approximation as shown in the figures above. 



Figure 6 

We see that even the “largest” approximation of Cd{z) is much smaller and cheaper 
to compute than Cd{z) and becomes more efficient as the degree of the polynomials 
increases. Since our main interest is in determining when some given z G C is a root 
of some polynomial in jCd, we observe that it would be significantly cheaper to check 
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if the “largest” approximation intersects some 5-ball about 0. 


In conclusion, by looking at the group structure of the set of Littlewood polynomi¬ 
als for some fixed degree d, and under the appropriate group operation, we may build 
easily computable approximations of its zero set. Computing the zero set of Cd in¬ 
volves finding roots of a very large number of high degree polynomials, which, while 
it is a straighforward procedure, it is computationally very expensive. Luckily, the zero 
set of Cd has enough symmetry to allow for a much cheaper approximation. 
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